Stability and dewetting of thin liquid films 
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Summary. The stability of thin liquid coatings is of fundamental interest in every- 
day life. Homogeneous and non-volatile liquid coatings may dewet either by heteroge- 
r~j— ^ neous nucleation, thermal nucleation, or spinodal dewetting. Wetting and dewetting 

Q | is explained on a fundamental level, including a discussion of relevant interactions. 

CZ2 , The chapter will also address the various dewetting scenarios and explain how the 

effective interface potential governs the behavior obtained for various stratified sub- 
strates and film thicknesses. 



^3 

£h ' 1 Introduction 

o ■ 

'To wet or not to wet' pQ is the question this chapter tries to answer. It con- 
cerns the fundamental aspects of the stability of thin liquid films or coatings. 
The liquid in question can be any homogeneous liquid^!, such as, e.g., an aque- 
ous solution, an oil, or a polymer melt. Whether or not a liquid wets a given 
' surface is not only of fundamental interest, but also of substantial techni- 

, cal importance. The authors have been approached by numerous companies 

with interesting problems concerning wettability, ranging from the automotive 
ly-^ , industry, the pharmaceutical and chemical industry, to the food processing, 

' printing, and textile industry. In principle, the general answer to all wettabil- 

ity problems is simple: a liquid wets a surface if it can gain energy by enlarging 
the interface with that substrate. To provide this answer for a specific system 
under study, however, may be quite cumbersome. Below we provide what may 
be seen as a map to reader to help with a first encounter with thin liquid film 
stability problems. For the sake of clarity, we restrict our considerations to 
' liquids on solid, homogeneous substrates. 



3 homogeneous on the length scale set by the thickness of the liquid film under 
consideration. In particular, demixing or segregation effects are thereby assumed 
to be negligible. 
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2 Experimental model systems for simple liquids 

In order to study the basic mechanisms of a certain class of phenomena, it 
is of central importance to identify suitable model systems. Since the main 
problems in the controlled preparation of thin films and surfaces lies in their 
inherent propensity to become contaminated by dust or other impurities, the 
first quantitative studies of surface forces in wetting films have been performed 
with cryogenic systems, in which impurities are naturally frozen out [21 EH]. 
However, it has been realized that these systems involve strong exchange of 
material between the liquid and vapor phase and thus are not well suited 
for studies of dynamic aspects, such as dewetting and structure formation 
phenomena. 

Thin laser-annealed metal films were candidates, [6] yet the same time, 
techniques of preparing systems which involve more complex liquids, such as 
polymer melts, in a well-defined and clean way became progressively available. 
As a result, the numerous benchmark studies of the basic mechanisms of 
dewetting have been performed with polymer melts (see e.g. the review articles 
of D.G. Bucknall[7] and P. Muller-BuschbaunrSi). Systematic experimental 
studies were begun somewhen around the 1990s QJjl [UJ H21 Q3] inspired 
by P.-G. de Gennes' theoretical work. [T5I RH] Polymer melts are on the one 
hand close to application (coatings, photo resist), yet on the other hand easily 
controllable in the experiments. Polymers such as polystyrene (PS) are very 
suitable model liquids since they have a very low vapor pressure in the melt, 
and mass conservation can safely be assumed. Moreover, they are chemically 
inert, non-polar, and their dynamics can be tailored by choosing different 
chain lengths and annealing temperatures. For molecular weights below the 
entanglement length (~17 kg/mol), the melt can safely be treated as simple 
(Newtonian) liquid for the low shear rates in dewetting experiments|16[ \T7\. 
Below the glass transition temperature T g , the films are glassy and can be 
stored for subsequent analysiqj. PS in atactic steric configuration moreover 
does not exhibit any tendency for crystallization. 

As substrates, typically hydrophobized Si wafers with their natural amor- 
phous Si dioxide layer (SiOx) have been widely used, the present chapter thus 
concentrates on Si wafers as substrates, which can be purchased with a very 
small surface roughness (rms roughness smaller than 0.2 nm). The wettabil- 
ity of these substrates can be tuned by chemically grafting monomolecular 
functional layers onto the native silicon oxide overlayer. For instance, hy- 
drophobization is often achieved by grafting a self-assembled monolayer of 
octadecyl-trichlorosilane 25, 26J onto the wafer. 

For dewetting experiments it is necessary to prepare a thin liquid film in 
a non-equilibrium state on a substrate. Usually, a thin polymer coating is 
prepared from a solvent solution by one of the following standard techniques: 



4 It is important to note that T g of thin films can be substantially different (mostly 
lower) than the bulk value. Q21 US I2TJ1 |2"T1 12"2"1|2"51 12"3] 
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spin coating, dip coating, or spraying. The solvent evaporates during the 
preparation procedure, leaving behind a smooth, glassy polymer layer. By spin 
coating a solution of polystyrene in toluene one can easily achieve a PS layer 
of thickness in the range of a few nanometers up to several micrometers. The 
roughness of the polystyrene layer is then similar to the underlying substrate. 
For experimental details and caveats of preparation, the reader is referred 
to the pertinent literature. [27l [13l EH HUJ [Ml HSJ [3D] The thickness of the 
films can be determined by standard techniques, such as ellipsometry[3i . The 
dewetting process can be studied by optical microscopy and/or atomic force 
microscopy (AFM). Most studies use AFM in non-contact mode operation to 
avoid any damage of the soft surface, which even allows to image the dewetting 
scenario in situ[32l [33], 

To induce dewetting, the films are heated above the glass transition tem- 
perature. Figure [1] shows a series of optical micrographs of an 80 nm thick 
PS film of molecular weight of 65 kg/mol ('PS(65k)') dewetting a silanized Si 
wafer. The series depicts a pattern formation process that is typical for most 
dewetting films, whether for a water film on a waxed surface, or a coating on 
a dusty substrate. 




Fig. 1. Pictures series taken by a light microscope: a 80 nm thick polystyrene film 
of 65 kg/mol molecular weight is dewetting at 135° C from a hydrophobized silicon 
substrate for approximately (from left to right) 2 min, 10 min, 30 min, and 100 min. 

The process of dewetting can be divided into three stages: in the early 
stage, holes are generated by a rupture process, c.f. Fig.fTk; in the intermediate 
stage, the radius of the holes increases, leading to hole coalescence, c.f. Fig.[T|D- 
c. In the intermediate stage, the focus is on the dynamics involved in the 
dewetting process, its impact on the hole profiles and on its influence on 
dewetting patterns. From the dynamics of hole growth [571 [3H 133 1351 [33 [551 
[551 [551 HO] as well as from the shape of the liquid rim surrounding the hole 
one can get information about the slip or no-slip boundary condition of the 
liquid close to the solid substrate [T71 2H [32] . In the late stage, the straight 

5 For all these techniques, the surface to be coated must be wettable by the polymer 
solution. Otherwise, the film must be prepared on another substrate, e.g. mica, 
and then transferred to the surface of interest. 
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ribbons that separate two coalescing holes decay into droplets due to the 
Rayleigh-Plateau instability 43J, c.f. Fig.[T]>d. Slight differences in the size of 
the droplets cause slight pressure differences, leading small droplets to shrink 
and large droplets to grow (Ostwald ripening[33]). Due to the small vapor 
pressure and the low mobility of the polymer molecules, the ripening process is 
extremely slow. Most experiments are stopped before the 'final state', i.e. one 
single drop on the surface, is reached. In what follows, we focus on the initial 
stage of dewetting: why does the polystyrene film dewet the hydrophobized 
Si wafer at all? The next section will be dedicated to the discussion of the 
various energy contributions involved, and their mutual balance. 



3 The energy balance 

In a Teflon© coated frying pan a stable oil layer can only be achieved if 
one pours enough oil into the pan. Gravity then stabilizes the thick oil film. 
If the oil film thins well below the so-called capillary length A CQp , capillary 
(intermolecular) forces dominate over gravitational forces and the film may 
dewet [IS]. The capillary length is given by X cap =y / iJi v / pg (for olive oil: X cap « 
1.7 mm, for water: A cap « 2.7 mm), where ai v is the liquid/vapor surface 
tension and p the density of the liquid. Therefore, dewetting of liquid films 
driven by intermolecular forces can only occur for 'thin' films, i.e. films thinner 
than A ca p. 

A liquid can be metastable and then needs a nucleus to induce dewetting. 
In the following we will focus on liquid films thinner than X cap and discuss un- 
der which conditions a thin liquid layer will be stable, metastable or unstable 
on top of a substrate. 



3.1 The Young equation 

Dewetting is a dynamic process that begins in a non-equilibrium situation 
[461 HZ1 HH1 H3 [501 Ell HH [52], namely the flat film on the surface and ends 
when reaching an equilibrium state, one droplet or a set of droplets 0. Thus, 
let us first start with a droplet, c.f. Fig.O 

A droplet on a homogeneous surface usually exhibits the form of a spherical 
cap and the tangent to the droplet at the three phase contact line includes an 
angle 6 with the substrate. In equilibrium, the angle is given by a balance of 
macroscopic forces: 



This is the Young equation [46] of 1805, where a sv and a s i are the solid/vapor 
and solid/liquid interfacial energies, and ai v is the liquid/vapor interfacial 
energy (or tension). 

6 Equilibrium situation actually is one single droplet, yet it normally takes too long 
to reach this state. Before, a network of droplets is formed. 
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Fig. 2. Sketch of a liquid drop atop of a solid substrate. Complete wetting is char- 
acterized by a contact angle 9 = 0, partial wetting by < 8 < n, and nonwetting 
by 6 = 7T (from left to right) 

For 6 = 0°, i.e. a sv ~ a s i > ai v the droplet will spread and will completely 
wet the substrate. For 0° < 9 < 180° one speaks of partial wetting, and for 
9 = 180° of non- wetting. All three cases are sketched in Fig. [21 In case of water 
(resp. oil), a surface is termed 'hydrophilic' ('oleophilic') if 0° < 9 < 90° and 
'hydrophobic' ('oleophobic') for 90° < 9 < 180° 0. In other words, the contact 
angle in Young's equation is determined by the free energies of interfaces 
between semi-infinite media. 

This is, however, not the full story. Although a first glance at Fig. [2] 
suggests that the free energy of a homogeneous film should be written as 
ffiim = °si + °iv , this approach neglects possible interactions of the two inter- 
faces (solid/liquid and liquid/vapor) with each other, across the liquid film. 
Such interactions may come up by virtue of van der Waals forces between the 
involved materials, which may span several tens of nanometers. In the vicinity 
of the three phase contact line, the interfaces can thus deviate significantly 
from a straight intersection at Young's angle (cf. Fig. [3]). We will come back 
to this point later, after having discussed in some detail the long-range forces 
of a stratified system. 

These forces are of similar importance for predicting the stability of a liq- 
uid film. Their dependence on film thickness determines not only the linear 
stability against small perturbations of the free surface, but also the scenario 
by which the initially uniform film is transformed into its equilibrium state, 
which consists of droplets of contact angle 9 on the surface [331 1SH 1SS] • For 
example, a droplet of photo resist on a semiconductor may exhibit an equilib- 
rium contact angle 9 — 20° and therefore wets the surface only partially. What 
does this mean for a photo resist film that was prepared by a non-equilibrium 
technique (spin coating) on top of the semiconductor? According to Young's 
equation, the film is not stable. How will the film decay into droplets and at 
what speed will this process take place at a given temperature? Will there be 
time to dry or cure the photo resist before it dewets? 

Depending on the thickness of the photo resist, the answers will be dif- 
ferent, though Young's contact angle 9 is always the same! Macroscopic and 



7 It should be noted that non-wetting has never been observed on flat surfaces, 
although some systems (like, e.g., mercury on glass) come quite close. 
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Fig. 3. Droplet profiles in the vicinity of the three phase contact line as expected 
for a) an unstable and b) a metastable situation. The profiles have been calculated 
for the case of a) a dsiOx = 191 nm on Si and b) a dsiOx = 1.7 nm SiOx layer on Si, 
the potentials of which are shown in Fig. [7] The dashed line marks the macroscopic 
contact angle 9 = 7.5° for a PS drop and the arrows mark the equilibrium wetting 
layer thickness h eq 

molecular terms describing stability conditions often were mixed up and have 
led to confusion. In what follows, we shall discuss which forces are to be con- 
sidered, and what typical length scales are involved. 

3.2 The effective interface potential 

Following the preceding section, we define the so-called effective interface po- 
tential <fi(h) = a s i +<7i v — <Tfii m . It comprises both short-range and long-range 
interactions and is defined as the excess free energy per unit area which is 
necessary to bring two interfaces from infinity to a certain distance. From its 
very definition, it is clear that 4>{K) — > as h tends to infinity. The excess free 
energy of an infinitely thick film is thus given by the sum of the free energies 
of its two interfaces, in accordance with intuition. 

For dielectric systems, there are only two relevant types of interactions, 
steric repulsion and van der Waals forceqj. 

(j)(h) = 4>(h)steric + 4>{h)vdW (2) 

In Fig. [4] three typical curves of 4>(h) are sketched to illustrate the general 
principle. Line (1) characterizes a film that is stable on the substrate, since 
energy would be necessary to thin the film. The equilibrium film thickness is 
infinite. The two other curves exhibit a global minimum of <f>(h) at h — h eq : 
Curve (2) characterizes a film that is unstable, whereas line (3) describes a 
metastable film. 



A detailed description for numerous different situations involving e.g. polar 
molecules or hydrogen-bonds or the interactions of a colloidal sphere interact- 
ing with a metallic surface can be found in the textbook of J. Israelachvili 48 . 
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Fig. 4. Sketch of the effective interface potential as a function of film thick- 
ness. Line (1) denotes the stable case, line (2) the unstable one and curve (3) the 
metastable case. The thickness of the stable wetting layer is termed h — h eq , and is 
typically in the order of some nm, the depth of the global minimum of <f> is named 
</>(h eq ). 



It is readily shown by a linear stability analysis [56] [5TJ [58] that if the 
second derivative of 4>{h) with respect to h is negative {(j)"{h) < 0), unsta- 
ble modes exist whose amplitudes grow exponentially according to exp(i/r), 
where r is the growth time that is characteristic for the respective mode. Fur- 
thermore, there is a characteristic wavelength X s of these modes the amplitude 
of which grows fastest and will therefore dominate the emerging dewetting 
pattern. 

The linear stability analysis also reveals that the spinodal wavelength X s 
is linked to the (second derivative of the) effective interface potential [55] [57] : 

The spinodal wavelength A s is the key to clearly identify a spinodal dewetting 
process in experimental systems^] 

One question is now to be answered: What links the effective interface po- 
tential to Young's contact angle? The effective interface potential can describe 
a non-equilibrium situation, yet Young's contact angle is only defined for the 
equilibrium. A. Frumkin has published in 1938 that [60] 

^2=co S 0-l. (4) 

9 The name 'spinodal dewetting' has been coined due to the analogy to spinodal 
decomposition of a blend of incompatible liquids, which occurs if the second 
derivative of the free energy with respect to the composition is negative [59]. As 
in dewetting, this mechanism is clearly distinct from heterogeneous nucleation, 
in which randomly distributed impurities dominate the emerging pattern. |13| 
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What is the consequence of Eq. and Eq. ^? Determining the spinodal 
wavelength A s as a function of film thickness h enables us to gain insight into 
the course of <f>" (h). By additionally measuring the equilibrium layer thickness 
h eq and the contact angle 8, it is possible to reconstruct the complete effective 
interface potential, [29l [301 E] c.f. Fig.[7l as will be shown later. 

The steric and the van der Waals part of the effective interface potential 
are characterized by different exponents and different interaction constants, 
therefore we will discuss the interactions separately: 

Steric repulsion and chemical interactions are relevant only within a few 
Angstroms of film thickness, and the resulting force is therefore termed 'short- 
range force'. The repulsion is due to overlapping electron shells and is typically 
described by a higher-order polynomial function and varies as l/h 12 , where 
h is the distance between the interacting bodies and is difficult to quantify. 
Considering two planar surfaces, this repulsion yields an interaction energy 
varying as 

C 

4>(h) steric = Tg j (5) 

where C is a constant characterizing the interaction strengtrF"). The reason 
for the lower exponent in the flat film geometry lies in fact that mutual in- 
teraction between all involved atoms or molecules have to be considered. For 
an extended derivation of the equations we like to refer to J. Israelachvili's 
approach^!. 

Van der Waals interactions characterize attractive intermolecular 
forces of quantum-electrodynamic origin. They arise from the variations of 
the zero-point energies of the collective electromagnetic modes of the system 
under study. The van der Waals energy between two molecules turns out to 
vary as l//i 6 in the non-retarded approximation. This approximation holds 
as long as the lateral dimensions of the system are much smaller than the 
wavelength of the electromagnetic fields at the dominant excitation energies. 
Considering retardation effects, the interaction falls off as l/h 7 . In the follow- 
ing, we use only non-retarded potentials since in the experimental systems, 
effects due to retardation are mostly camouflaged by experimental error bars. 
Nevertheless, in some special cases including cryogenic systems, retardation 
effects have been reported [3] and a recent theoretical study [H] discusses their 
possible influence on thin films in the context of the full theory of Dzyaloshin- 
skii, Lifshitz, and Pitaevskii [53] , 

Considering again two planar surfaces, the non-retarded interaction yields 

= -i2^a> ( 6 ) 



10 This definition of tf>(h) s teric here is such that C is positive and thus the interaction 
always repulsive. 

1 Here, chapter 10 in the textbook of J. Israelachvili [48] is very helpful. 



Stability and dewetting of thin liquid films 9 

where A is the Hamaker constant [64] and (f>(h) is the energy per unit area. A 
retarded interaction is described as <j>{h) ~ A/h 3 . 

The van der Waals forces of two media interacting through vacuum are 
always attractive, and A is positive. If the vacuum yet is replaced by a third 
medium (e.g. a liquid film), things get more involved: In a system consisting 
of three media 7711/7713/7712, it may happen that mi attracts 777.2 stronger than 
777,3 attracts 777,2. Thus 7712 is 'repelled' by 7773, this means that the Hamaker 
constant is negative. To describe the entire system, the Lifshitz theory[65] has 
to be applied. Here, the interacting bodies are treated as continuous media 
and the atomic structure is ignored [65l [48l [62] . Rather, bulk properties as 
the dielectric constants and the refractive indices are used to calculate the 
Hamaker constant. 

For the sake of completeness, we cite here the formula @ for the Hamaker 
constants from the book of J. Israelachvili [4"8] . which has proven to be very 
useful. It is valid for two media 1 and 2, interacting through media 3. All media 
are taken as being dielectric with a single electronic absorption frequency v e , 
which is typically in the range of 3 • 10 15 Hz; 77 is the refractive index of the 
medium in the visible and e the dielectric constant (i.e. n 2 = £j, taken in the 
visible spectral range) 

A = A u =q + A u>0 

4 \ £ 1 + e 3 J \ £ 2 + £3/ 

ihv e (nl — n§)(n| — 773) 

+ 8V2 y/(nl + WTWl + WWWx + n%) + V(n2+n§)} 



(7) 

(8) 



The Hamaker constant A will then give the strength of the van der Waals 
forces between the two interfaces solid/liquid and liquid/air. 

For air as medium 1, Si as medium 2 and polystyrene as medium 3, Eq. [5] 
give0 Asi = —2.2(5) • 10~ 19 J, replacing Si as medium 2 by silicon dioxide, 
the Hamaker constant is Asiox — 1-8(4) • 10~ 20 J. For a silane (OTS) covered 
Si wafer, the Hamaker constant for the interaction with PS is calculated to 
be Aots = 1.9(3) • 10~ 20 J. For A < (A > 0), the system can gain energy 
by enlarging (reducing) the distance h between the surfaces, in other words, 
a polystyrene film of thickness h is stable on Si, since A < 0. 

In Fig. the dotted line represents the van der Waals potential, 4>(h) v dw, 
as given by Eq. [6] for a polystyrene film on a Si wafer, where A$i is negative. 
The potential therefore is positive, purely repulsive £3 and the PS film will be 

12 Confer chapter 11, Eq. 11.13 of Ref.gg] 

13 The error of the Hamaker constant is a result of the uncertainties of the optical 
properties of the involved materials as found in literature. The error for the last 
digit is given in brackets. 

14 The terms 'attractive' and 'repulsive' can be misleading if thinking of the thin 
film, yet the terms are chosen for a system of two media (in our case solid and 
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stable. However, (f>(h) vc tw f° r PS on an infinitely thick SiOx layer (solid line) 
is always negative, purely attractive and the PS film will be unstable, since 
the Hamaker constant AsiOx is positive. 




polystyrene film thickness [nm] 

Fig. 5. Long-range part of the effective interface potential <f>(h) as function of PS 
film thickness h for different SiOx layer thickness ranging from nm (dotted line) 
to infinity (solid line), calculated with the formula given in Eq. (|10|l . The Hamaker 
constants were calculated from the optical constants of the involved materials |48|. 

Due to the 1/h 2 dependence of the potential, the van der Waals forces are 
long-range forces and act significantly on distances up to about 100 nm. In 
stratified systems with more than one layer between the two half spaces of 
media 1 and 2, all mutual interactions have to be considered. The exact calcu- 
lation of the van der Waals potential may therefore be quite cumbersome. This 
is an important point, since there is a number of publications that disregard 
this aspect and claim that theory and experiment do not match. A system like 
a PS film in air, spun on a thin silicon oxide layer on top of a Si wafer involves 
two thin layers (three interfaces, abbreviated as air/PS/SiOx/Si) and can be 
described by two Hamaker constants. So far, no discrepancies between the- 
ory and experiment have ever been observed when the formulas given above 
had been used. One on the most impressive successes is the quite accurate 
prediction of the critical wetting temperature of pentane on water [66j . 

Assuming additivity of forces, the system air/PS/SiOx/Si can be charac- 
terized by a summation of the van der Waals contributions of each of the 
single layers with thickness h of the PS and dsiOx of the SiOx layer: 



air) interacting through a third one, here the liquid layer. Hence attractive and 
repulsive are meant for air being attracted to the solid surface or 'repelled'. In 
case it is repelled, the effective interface potential is repulsive for the respective 
thickness of the liquid layer and the liquid wets the solid. 
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, /,s A Sl0x A Sl0x — A Si 

<t>vdW{h) = - + , , (10) 

Yl-nh 1 I2ir[h + dsiOxY 

where AsiOx and Asi are the Hamaker constants of the respective system 
air/PS/SiOx and air/PS/Si. Although not exact^ the concept allows to cal- 
culate the interaction energies even for stratified systems. With the help of 
Eqs. © and (fTUl) , the van der Waals potential of the experimental system is 
accessible, if the SiOx layer thickness is known and the Hamaker constants 
are calculated as shown before. Fig. [5] depicts the van der Waals potential as 
gained from Eq. (fTU)) . Clearly, the potential is influenced by the thickness of 
the silicon dioxide layer. 

So the widely used system air/PS/OTS/SiOx/Si can be described by 

rh (h\ - AqTS AqtS — AsiOx AsiQx — Asi , . 

<p vd wW Unh2 + 12<h + doTs)2 + + ^ + , (llj 

yet the long-range potential 4> v dw is n °t much different to a system with an 
SiOx layer of thickness dors + dsiOx , since the optical properties of the OTS 
and the SiOx layer are very similar. The difference in wetting behavior will 
only be obvious if the short-range potential is added. 

To conclude so far, <p(h) combines short- and long-range interactions. 
Short-range interactions are difficult to quantify, yet van der Waals inter- 
actions can be captured by the optical properties of the involved materials, 
even for stratified systems, where the suitable van der Waals potential must 
be taken. If <j>"{h) < 0, spinodal dewetting is possible and a characteristic 
wavelength A s should be imprinted on the dewetting pattern. 



4 Experiments: Linking the effective interface potential 
to macroscopic properties 

In this section experimental results shall be compared to the predictions made 
by the effective interface potential. Typical experimental dewetting patterns 
in the system air/PS/SiOx/Si are shown in Fig. [6l What can we learn from 
analyzing the dewetting films? I) By AFM, the contact angle 9 of droplets in 
the late stage of dewetting can be measured. II) The equilibrium film thickness 
h eq can be determined, e.g. by ellipsometry. Ill) The spinodal wavelength A 

15 It should noted here that numerous studies have tried to use only one so-called 
'effective Hamaker constant' for stratified systems using combining rules of the 
form A132 ~ (V A22 — y/Aix) ( VA22 — V ^33) • However, these rules can only achieve 
good results under severe limitations, e.g., if the zero-frequency contribution A„ = o 
in Eq. (JSJ is negligible 48, . Strictly speaking, they are plain wrong, and should be 
used, if at all, with utmost care. It is not surprising that many of the studies which 
have used them could not quantitatively reconcile the theoretical description with 
the experimental results. 
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and IV) the growth time r of the spinodal pattern can be determined. V) The 
morphology of the pattern. VI) The growth of the size of the holes and VII) 
the form of the liquid front as function of time. 




spinodal dewetting thermal nucleation heterog. nucleation 



Fig. 6. a)-d) AFM images of dewetting PS (2k) films on Si wafers with variable Si 
dioxide layer and PS film thickness. Scale bars indicate 5 /im, z-scale ranges from 
(black) to 20 nm (white): a) 3.9 nm PS film on a Si wafer with dstox =191 nm. 
The inset shows a Fourier transform of the image, b) 3.9(2) nm PS, c) 4.1 nm PS 
and d) 6.6 nm PS on Si wafers with dsiOx =2.4 nm. (The statistical analysis of the 
distribution of hole sites in cases (b) to (d) was performed on larger sample areas.) 

Dewetting patterns like the one in Fig. [6^i clearly exhibit the spinodal 
wavelength. Fig. [7^ comprises experiments with different PS film thicknesses 
on Si wafers with dsiOx — 2.4 nm (open symbols) and with dsiOx — 191 
nm (filled symbols). From the data of X(h) data points for <fr"(h) can be 
gained. Now, fitting the second derivative of Eq. (flTj)) to the data, the Hamaker 
constants and the short-range interaction constants can be obtained as fit 
parameters and the entire effective interface potential can be inferred. It is 
plotted in Fig. [7b. 

The best fits [25] |3D] are achieved for A SiJU = -1.3(6) • 1(T 19 J and 
AsiOxjit = 2.2(4)-10~ 20 J. The values for A nicely match the values calculated 
from optical properties of the media involved, as described by Eq. Ejl For C 
we find CsiOxjit = 6.3(1) • lO" 76 Jm 6 and CqtSJU = 2.1(1) • 1CT 81 Jm 6 . In 
the cited references, a detailed description of the fitting and reconstruction 
procedure can be found. 

The second derivative <j>"{h) is plotted in in Fig. [8^ for three SiOx layer 
thicknesses. An arrow marks the zero in the respective second derivative, cor- 
responding to the relevant inflection points of the curves displayed in Fig. [JJd. 
Spinodal dewetting is possible only for films thinner than indicated by the 
arrows: PS films on Si wafers with an 1.7 nm (2.4 nm) thick oxide layer are 
unstable for PS film thicknesses below 3 nm (4 nm). For larger thicknesses, 
the films are metastable and need to overcome a potential barrier in order to 
finally reach the global minimum of <f>. The maximum in (p, which is visible 
in Fig. [JJd for the broken curves, represents only part of this barrier, since a 
nucleus for dewetting as it is formed, e.g., by thermal activation, is a localized 




Fig. 7. a) Spinodal wavelength as function of PS film thickness on Si wafers with 
dsiOx = 2.4 nm (open symbols) and with dsiOx = 191 nm (filled symbols), b) The 
effective interface potential 4>{h) for three different SiOx layer thicknesses. 



structure and thus involves excess surface energy as well [55]. It is readily 
checked that the energy required to form such nucleus is in almost all cases 
large as compared to kT, such that thermal nucleation plays no role in sys- 
tems of practical interest. This can be illustrated nicely with the PS films, 
as shown in Fig. [Hb- PS films on thick 191 nm SiOx (solid line in Fig. [7^) 
are unstable for all relevant film thicknesses. On 2.4 nm of oxide, however, 
a sign reversal is observed at a film thickness of about 4 nm (white arrow 
in Fig. [8^,). Only very close to this point, the height of the potential bar- 
rier vanishes, and homogeneous nucleation by thermal activation is possible 
[68l [45l [69] . However, nucleation can as well, and usually does, proceed by 
means of localized defects in the film. A defect, be it in the molecular texture 
of a polymeric film or just a small dust particle, may remove the potential 
barrier locally and thus induce dewetting. This rupture mechanism is termed 
'heterogeneous nucleation' [59l|67], in analogy to defect-mediated nucleation 
in demixing scenarios. 

The fact that the effective interface potential can be directly inferred from 
macroscopic quantities, such as the spinodal wavelength, is one reason why 
experimentalists are seeking for spinodally dewetting regimes in thin liquid 
film systems. However, if a company is asking why their coating does not stay 
stable on a surface, it is convenient to rather have a stability diagram which 
allows to look up where the system will be in a metastable or an unstable state. 
A stability diagram for the system air/PS/SiOx/Si is shown in Fig. [SJa. The 
solid lines of the diagram base on Eq.[5]with a long-range potential as in Eq.fTUl 
and separate the spinodal dewetting (unstable) regime, where <j>"(h) < 0, 
from in the regime of heterogeneous nucleation, characterized by <f>"(h) > 
(metastable regime). Thermal nucleation is possible for <fi"(h) = 0, this 
line separates the two regimes. Experiments that exhibit spinodal dewetting 
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patterns are indicated by open symbols, those of heterogeneous nucleation 
(randomly distributed holes) are indicated by solid symbols^]. A star marks 
the set of parameters where thermal nucleation was observed. Note that in 
the unstable regime, heterogeneous nucleation from localized defects is also 
possible and indeed is sometimes observed. In fact, it usually dominates if 
the experiments are not performed in extremely clean conditions. Spinodal 
dewetting, however, can only take place in the unstable regime. 




Fig. 8. a) Second derivative of the effective interface potential as function of 
film thickness for three different SiOx layer thicknesses as fitted to the data of 
Fig. [7] b) Stability diagram of PS films on top of Si wafers with variable oxide layer 
thickness. 1291 



We can state here that the effective interface potential indeed is the key 
for wetting properties of liquids on surfaces. It can easily be obtained via 
the optical properties of the involved media and by taking into account the 
adapted van der Waals potential for stratified systems. It is worth noting that 
metastable potentials can only be present if there is more than one Hamaker 
constant involved, one with negative and one with positive sign. 

Coming back to the company - a stability diagram for their specific system 
will reveal the way to a solution of their wettability problem: If their coating 
system is located in the metastable regime (in the example shown, it depends 
on the silicon dioxide thickness and the polystyrene thickness), the advice 
is to reduce the number of possible nucleation centers, e.g. by creating a 
cleaner environment. If the coating system yet is in the unstable regime, the 
introduction of an adhesion promoter is necessary. 



As explained above, thermal nucleation can be only observed if the nucleation 
barrier is of order kT. On the scale of Fig. [8]d, this is fulfilled only in a region 
smaller than the width of the line denoted 'thermal nucleation'. Alleged observa- 
tions of thermal nucleation in a wider range [70] appear very questionable. 
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As mentioned at the beginning, the shape of a droplet close to the three 
phase region is influenced by the long-range intermolecular interactions. By 
AFM, this region can be characterized. For a detailed discussion we like to refer 
to the work of T. Pompe et al. [71] l72 ] . In Fig. [3J two profiles are sketched that 
were calculated based on the knowledge of the effective interface potential. 
Fig. EH is a characteristic curve for an unstable potential as shown for dsiOx — 
191 nm in Fig. [Jji, whereas a profile like Fig. [3Jd is expected for a metastable 
potential like for dsiOx = 1.7 nm. The experimental AFM results (not shown 
here) corroborate the theoretical expectations [75] . 

One question is open up to now: are spinodal dewetting patterns always 
as easy to detect as in Fig. [6^? The answer is 'no': patterns which at first 
glance appear like the one shown in Fig. [6Jd can be generated due to a spinodal 
process as well |6] . The next section we will briefly explain how to distinguish a 
dewetting pattern of a system in the metastable state from one in the unstable 
regime and also address thermal nucleation. 



5 Characterizing experimental dewetting patterns 

The three rupture scenarios spinodal dewetting, homogeneous nucle- 
ation, and heterogeneous nucleation give rise to specific dewetting pat- 
terns. Vice versa, characterizing the experimental dewetting pattern can help 
to identify the rupture mechanism and to infer the effective interface poten- 
tial. Fig. [TJ and Fig. [5] show experimental examples of films dewetting via 
different rupture mechanisms. Theoretically, the distinction between nucle- 
ation and spinodal dewetting appears quite clear: Vrij |56j proposed already 
in 1966 that a spinodal rupture of a free liquid film results in a dewetting 
pattern of 'hills and gullies' with a preferred distance A s after a certain time 
of rupture r. 

Experimentally, the rupture time r is difficult to measure since the hole 
which forms as the film has thinned to zero thickness must have a certain 
size to be observable. Experimentalists thus concentrated instead on the ev- 
idence of a preferred wavelength \ s observable in their systems^ If, how- 
ever, the holes are randomly (Poisson) distributed, they are assumed to stem 

17 Note that a preferred hole distance \ s has to be distinguished from a mean hole 
distance: In some studies it was found that the mean hole distance scaled with the 
film thickness as expected for spinodal dewetting [101174] . Therefore, the scenario 
shown in Fig. [Ik was often regarded as a typical spinodal dewetting scenario. 
However, it turned out later that this was not correct: [T3] first of all, the system 
is metastable for the film thicknesses studied, as an analysis similar to Eq. [TO] 
would have shown. Secondly, the typical time scale for spinodal dewetting docs 
not fit the theoretical expectation (rupture time r k /i 5 ). 
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from heterogeneous nucleation, reflecting the fact that defects typically 
exhibit random statistics. Although it was generally accepted later that (het- 
erogeneous) nucleation from localized defects is the reason for the dewetting 
scenario shown in Fig. []J,, the very nature of the nucleation defect mostly re- 
mains unclear: it might be dust particles or any other chemical or physical in- 
homogeneities. In some holes, by light microscopy or atomic force microscopy 
(AFM), a nanoscopic object could be observed right in the center of each hole. 
Assuming the object to be a dust particle, the first trial was to reduce the 
number of nucleation centers by improving the preparation conditions. The 
number of holes, however, could not be reduced below a certain level, which 
suggests that the physics of hole nucleation in polymer films may be deeper 
than a mere effect of 'dirt'. It has been shown experimentally that stress inside 
the thin films (stretched entanglements [75] ) caused by the preparation of the 
film out of solution plays a significant role and can cause holes. Details can 
be found in the studies of Seemann et al. [73J and Reiter et al. [3J5] . 

The experimental distinction between a hole pattern with preferred hole 
distance and a pattern with randomly distributed hole sites is far from being 
obvious. The thicker the films are, the weaker is the driving force, and the 
longer is the growth time r of the spinodal mode[56j[57l[30] (typically, r oc h 5 ) 
and can easily exceed experimental time scales. For thicker films, dewetting 
by heterogeneous nucleation may therefore be quicker and can suppress a 
spinodal pattern [76] [77]. Moreover, chemical heterogeneities can locally cause 
a change in <fr and therefore the rupture conditions of the sample may vary 
from spot to spot leading to a less-ordered dewetting pattern. This effect is 
more pronounced in thicker films due to the small driving forces and the large 
growth time r. Hence a two-point correlation function (a radial pair correlation 
function g(r) or a Fourier transform) might not be sensitive enough to detect 
the correlations between the hole sites, or the statistics is too poor. Then, more 
powerful tools have to be applied. Minkowski functionals - based on integral 
geometrical methods - have shown to be a versatile method to track down 
higher order correlations. There is no room here to dwell on this method 
in any reasonable detail. The reader is referred to the pertinent literature 

HI H3J [FS [31 [33]. 

Let us come back to thermal nucleation once more. As discussed above, the 
potential barrier for nucleation of a dry spot may in principle be overcome by 
thermal activation, provided the system is sufficiently close to the sign reversal 
of 4>"{h). The characteristic feature of this scenario is a continuous breakup 
of more and more holes, whereas nucleation from defects causes holes that 
emerge only within a sharp time window [28] . Fig. [Bp depicts an example for 
a 4.1 nm thick PS films on a wafer with a 2.4 nm SiOx layer: Holes of different 
sizes are observable. For that system, the stability diagram of Fig. 0b reads 
that </>"(4.1 nm) ss 0, which is another strong indication that the theoretical 
predictions corroborate the experimental observations and are able to capture 
the wettability of dielectric systems. 



Stability and dewetting of thin liquid films 



17 



a) b) c) d) 




i 1 

8 jJITl 



Fig. 9. Dewetting morphology of a 3.9(2) nm PS(2k) film on a Si wafer with a 
191 nm thick SiOx layer as recorded by in situ AFM. Up to about 5000 s (e) the 
temperature was held constant at 53° and annealing times are given in the pictures. 
Afterwards, the temperature was successively increased to 100°. The scan size in 
(a)(d) is 1.5 m. In (e)(g) larger scans were made to check possible damage of the 
sample by the AFM tip. Scan (e) was taken some minutes before scan (d), whereas 
scan (f) was recorded some minutes after scan (d). Scan (g) characterizes the end 
of the dewetting process and Ostwald ripening of the droplets slowly sets in. 

6 Dynamics of spinodal dewetting 

At the beginning of spinodal dewetting studies, only one system was known to 
dewet spinodally: thin gold films on top of quartz substrates jSJ- Metal films, 
however, are not so easy to deal with as compared to polymer films, since 
the films need to be annealed by a laser. Hence, the time scale is in the ns 
range, crystallization plays a role and the Hamaker constants are not so easy 
to calculate. Polymer films such as polystyrene (PS), however, are dielectric 
and Hamaker constants can easily be determined via Eq. [5J 

In situ AFM studies on a dewetting PS film on a Si wafer with a thick 
SiOx layer reveal the dynamics of the structure formation process, c.f. Fig. |H1 
Fig. [10] depicts results of Fourier transforms of the AFM scans. Clearly, one 
mode is amplified fastest. The amplitude of the growing unstable mode, shown 
in the inset, grows exponentially with time, as expected from theory. Further 
investigations show that at the onset of the dewetting, no mode is specifically 
selected, rather, the amplitude reflects the roughness of the film. Comparing 
the results to deterministic simulations. [32] the pattern formation process and 
the morphology of the dewetting structures match very well, yet the tempo- 
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ral evolution of the morphology slightly differs. It turned out that thermal 
fluctuations accelerate the dewetting dynamics in the experiments. |69[ 133] 
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Fig. 10. Results of Fourier transforms giving the power spectral density of in situ 
AFM pictures at T — 53° C, on of which is shown in Fig. [Bk. Note that the time 
intervals between the curves are not constant. The inset depicts the amplitude of 
the undulation as a function of annealing time. The first data point at t < s gives 
the roughness of the PS film surface at room temperature (RT) as revealed from a 
Fourier transform. The solid line is a fit of an exponential growth to the data, as 
expected from theory 57 . 



7 Conclusion 

The wettability of a substrate is a delicate interplay of forces. For dielectric 
systems, the wettability can be described successfully by the effective interface 
potential </>, which is a sum of short- and long-range interactions. The latter 
are dominated by van der Waals interactions. Their strength can be obtained 
by calculating the Hamaker constant via the optical properties of the involved 
media. For stratified systems, the additivity of forces can be assumed. 

Rupture mechanism, pattern formation, morphology and dynamics of 
dewetting are all governed by the effective interface potential and the ex- 
periments can corroborate the theoretical expectations. 

The knowledge of the effective interface potential therefore allows for a tai- 
loring of the wettability on demand and is a great tool for further fundamental 
or applied studies. 
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